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ABSTRACT 



This thesis is concerned with the steady state response to a step 
input of a non-linear servomechanism having viscous friction and a load 
position unity feedback loop. A gear train with backlash between re- 
silient gear teeth was located between the motor and load. 

The physical equations describing the system are presented and 
adapted to a form for the Control Data Corporation l6o4 digital com- 
puter. The computer program used in this analysis is presented along 
with its flow diagram and description of its operation. 

The results of this thesis are presented in several forms using 
the parameters of system damping coefficient, distribution of friction 
and inertia, backlash angle and coefficient of restitution. The re- 
sults are applied to sample problems in illustration of their use for 
design and analysis. 

The authors wish to express their appreciation to Dr. George J. 
Thaler of the Department of Electrical Engineering, and to Dr. William 
Wainwright and Mr. Edward N. Ward of the Department of Mathematics for 
their assistance and encouragement in completing this work. 
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Tm 
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J 

Pm 

ft 

A 
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(ONE)* 

(THETAL) 

(THErADL) 

(kmotconst) 



(RHO) 

(MOINERT) 

(JLOADVAL) 

(FMOTOR) 

(fload) 

(delta) 



Angular step input (radians). 

Desired load displacement at steady 
state 

Output displacement of load (radians). 

Output velocity of load ( radians') 

' sec ' 

Error in radians. E- ~ 0 C 

Constant of proportionality relating 
position error to torque developed by 
the motor. 

Developed torque of motor. 

Torque of load. 

Torque of load referred to motor. 

Number of teeth on gear 1 
Number of teeth on gear 2 

Inertia of motor measured at motor. 

Inertia of load measured at load. 

Total inertial measured at motor 
J = Jm + 0\_ 

Friction of motor measured at motor. 

Friction of load measured at load. 

Backlash of gearbox measured in radians 
at the output. 

Slope of phase trajectory (combined 
system). $ cU c = £ 

-5 d ^ 

Slope of phase trajectory (load floating 

free) - A/J =_£ 

Velocities of motor and load after impact. 



*( ) Terms indicate computer mnemonics referred to in Sec. 3, Computer 
Program Development. 
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System damping coefficient. 
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System natural frequency. 




(OMEGANSQ) 


System natural frequency, squared. 


f*. 


(FLFTPRIW) 


Load friction ratio with respect 
to friction of the system. 
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(inertrat) 


Inertia ratio. 


^ L_ 


(thetam) 


Displacement of motor (radians). 


©vr\ 


(thetaem) 


Velocity of motor /radians\ 

\ sec ; 



1. Introduction and Background 

The problem under investigation was the response to a unit step 
input of a nonlinear servomechanism having viscous friction and load- 
position unity feedback loop. A gear train with backlash between the 
resilient gear teeth was located between the motor and load. A de- 
termination of the possible existence of a limit cycle with variation 
of several parameters was made. The variable parameters considered 
were motor , load, and system time constants, backlash and the co- 
efficient of restitution of the gear teeth. If a limit cycle existed, 
a further investigation wad made into the change of the size of the 
limit cycle with variation of the above parameters. 

If a gear train is required between the motor and load in a 
servomechanism, the gear train may be treated as ideal with perfect 
meshing of the gears. In this case the components of the system are 
joined at all times and the system may be described by a single linear 
differential equation. The classical linear solution in either the 
frequency or time domain is the result. 

Present production methods cannot meet the requirements of an ideal 
gear train, having perfect engagement of the teeth, nor is the ideal 
gear train desirable from the standpoint of wear on the gear faces. 

The result is that the practical gear train has separation of the gear 
teeth when the velocities of the motor and load are different. The re- 
sponse of a servomechanism with backlash is in effect discontinuous or 
nonlinear. The net system acts as three individual but dependent systems 
with the boundary conditions of each being specified by certain physical 
laws . 

Two accepted methods of analysis are available for solution of the 
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problem. The first is the solution by describing function methods 
while the second is by analytical methods. Describing function 
methods are primarily frequency response techniques and depend upon 
the assumption that the nonlinear element may be considered as linear 
over the range of consideration. A further assumption is made that 
the input to the nonlinear element is a pure sinusoid with no harmonic 
frequencies. The describing function method is primarily concerned 
with steady state results. 

The analytical solution, which is the second method of considering 
the problem is primarily a ^transient response technique concerned with 
solution in the time domain. This approach also requires certain 
assumptions concerning the behavior of the system. The assumptions 
used in several investigations are pointed out in the following para- 
graphs, while the assumptions of this investigation are outlined in 
Section 3> Development of Equations. The validity of the analysis by 
any method depends on the accuracy with which the real system has been 
defined. 

Chestnut and Mayer, Ref. a, have considered the backlash problem 
using describing functions and supported the findings by analog com- 
puter studies. The main difference between this cited work and previous 
describing function analyses was the simulation of backlash by a dead 
zone between the motor and load and springiness in the interconnecting 
shaft . 

Various analytical approaches have been made to the problem of 
backlash, Lutkenhouse, Ref. b, used graphical methods in the phase 
plane to solve several cases of plastic impact between the gear teeth 
of a second order servomechanism. Pastel and Thaler, Ref. c, developed 
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analytical equations for the existence of limit cycles using phase 
plane equations for the case of plastic impact and no load inertia. 

Knoll and Narud, Ref. d, investigated limit cycles in the phase plane 
using an analog computer. The investigation of Knoll and Narud 
covered a wide range of parameters for the case of plastic contact 
between the gear teeth. 

To analytically describe the total system response of a second 
order servo with backlash, three differential equations and one or 
more algebraic equations are required. The system as a whole may be 
treated as being piece wise linear. One differential equation is used 
to describe the entire system when the gear teeth are in contact and 
the motor is driving or braking. When the gear teeth are not in con- 
tact, two more differential equations are required, one for the motor 
alone and one for the load alone. The boundary conditions of the dif- 
ferential equations are determined from the solution of one or more 
algebraic equations, expressing the laws of conservation of momentum 
and energy. For plastic impact with no bounce of the gear teeth, the 
law of conservation of momentum must be satisfied. For perfect elastic 
contact between the gear teeth the law of conservation of momentum and 
the law of conservation of energy must be satisfied simultaniously. 
Intermediate cases between perfect elastic and perfect plastic contact 
can fulfill only the law of conservation of momentum. However, total 
energy accounting may be made in the intermediate and perfect plastic 
impact cases. 

New, Ref. e, adapted the differential equations of system motion 
and the law of conservation of momentum (plastic contact) to a digital 
computer analysis of a second order servomechanism with backlash. The 
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digital computer was chosen as the method of solution for this investi- 
gation primarily "because of the flexibility of data presentation and the 
ease with which cyclic parameter variation could be obtained. 

In later sections of this thesis, the differential equations for 
the system and the algebraic equations for the impact boundary conditions 
are developed. The computer program and its flow diagram and modes of 
operation are pointed out. The results of the solutions of problems 
are presented in several forms and observations are made in Section 
Results and Discussion. Sample applications of the results to mechani- 
cal systems are presented fn Section 6, Application of Results. 

At this writing, associated work in the area of transient response, 
primarily peak o-. ershoot and settling time, of a second order servo- 
mechanism with backlash and plastic and elastic impact is being pre- 
pared as a thesis by C. E. Andrews and R. A. Kelley at the U. S. Naval 
Postgraduate School. 
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2. Development of Equations 

The equations developed were those of a phase plane analysis. 

Prior to defining the net system equations, the assumptions under which 
the analysis was made will be stated. When necessary these assumptions 
will be amplified and referred to later in this work. 

It was assumed that: 

1. The gear teeth were initially in contact and the initial 
conditions of the system were all equal to zero. This was 
later proved to be an unnecessary limitation for the study 
of steady state response. 

2. Plastic deformation of the gear teeth during steady contact 
and impact and any torsional deformations of driving shafts 
are negligible. 

5. The inertias of the gears and drive shafts are considered 
as part of the load or motor inertia depending on their 
attachment in the system. 

4. The law of conservation of energy was completely satisfied 
in only the perfect elastic case by maintaining the total 
mechanical rotational energy of the system constant at the 
instant prior to and after impact. The law of conservation 

of momentum is satisfied in plastic, elastic, and intermediate 
cases. When the law of conservation of momentum is the only 
equation required to be satisfied, the energy lost from the 
system is dissipated in the heat of infinitesimal deformations 
of the gear teeth. 

5. The gear teeth are in contact only instantaneously during 
impact for the elastic and intermediate cases (excluding 
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plastic) and that the impulse torques of drive, friction and 
hearing supports, etc., are zero during impact, e.g. 



6. The coefficient of restitution of the two opposing gear 
teeth is the same or is described by an equivalent co- 
efficient if the two gear teeth are of unequal coefficients. 

7. Backlash is assumed to be equal at all points on the gear 
circumference. Backlash is measured at the output shaft. 

A block diagram of the' system considered is presented below. 
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The equations for the motor and load in contact are: 
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The equation for load alone with gears not in contact is: 

J L_ Oc + £l. ^C. ~ C) 

OV 0C 4~ 9|_ 0c - O 
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The equation for the motor alone with the gears not in contact is. 

Tw Qyv\ 4* f yy) ~ 



The method of analysis used by New, Ref. e, is used in order to 
be able to examine the system independently of its Co n 
Defining first RJn t = t* 



and differentiating 
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dt 1 - (dt*ndtj (dt*) 1 

d 2 &c - £h* and £L~ h) n 

cdt^r 

tnnalUl Oc.= &,t(x) A - Sc di* 

1 Cfc-tr 

^ cli = C dc dt* or 6t = &t 

Using the equations of the combined system and the load alone, 

6 C + Z ?0 l)/^ &C. + 6)yl 2 = d3v^ 

s S s> t e w> 

' * 4 

6c, t it &C. ~ & load 
Tl. *>***■ 

and making the indicated substitution to a transformed (*) coordinate 



system for the system equation 



Ud'A 6* + t2^(jOyo Sc. +■ U)y\ 6c = dj n &f>„ 

Sc ^§6 1 + at = e* 



and introducting the slope equations of the phase plane 

ft/s S t &C, = <9 results . 

If the load equation is first put in the form for the phase plane 

a/j_ t £±_ &c, ~ o 

0 i— ' / 

and then transformed to the (*) coordinate system where q Q 



ft/ L ^ f u ^ /> = & 



results . 
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The same result may be obtained by setting 4), =1. The results of such 
a transformation require inverse scaling for practical application, 
examples of which given in Section 6. It is noted that the slope of 
the load-free equation is mathematically the same whether in the phase 
plane or in the transformed phase plane. It is pointed out at this 
time that the equations later developed to satisfay the laws of con- 
servation of momentum and energy are independent of the system natural 
frequency. 

To establish equations for the law of the conservation of momentum 
and a relationship satisfying the law of the conservation of energy, 
the following schematic is used. 




It is seen that ^ ^ \ ) & C. ~~ ^ I 

_/ , 

and since P - 1 an equivalent schematic results 




By representing the inertia of the load and motor as inertia of the 
gears the following figure is obtained at the instant of impact. 
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The torque equation may be expressed \ = T9 - Tdj> 

(X t 

Impulse is equal to the time rate of change of momentum. The ex- 
pressions for the rates of change of momentum of the gears treated 
separately may he written: 
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where the primed velocities are those following impact. 

If the momentum is to he conserved in the system the impulse 
functions are equal and opposite. 

3J2L ( M ( K - & <) 

I * 

and 3m J_ ( - 6^) - ( 6- c - ) 

P 

from substitution of the equations 

r .' = r» , ©c = '(° 

The assumption that all other impulse functions are zero at the instant 
of impact is restated at this point. 

A definition of the coefficient of restitution will he made for 
purposes of this work: 

•e. - — ( ■=. — ( P 6^ - ^e, ) 

It will he shown that this definition of e will satisfy the law 
of conservation of energy. Consider the case of e=0, f 

6--( pi'- 

in which case for plastic impact (e=0) the gears are moving at the 
same velocity following impact. 

The expression for the coefficient of restitution is examined for 
the case of perfect elastic contact, e=l, for which the law of conser- 
vation of energy is satisfied. 

Solution of the restitution equation for e=l yields 
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and factoring 




The momentum equation which is independent of e is ; * ^ 

lift, f -- (6' c - 

p ^ u V - » / v / ' / * \ 

or rearranging ^ J “ ( £>C~ j 

The factored energy equation is now divided by the momentum 

equation to yield ^ ^ ~ / "& £ ~h &£, ) 

which checks with the definition of the coefficient of restitution 

for e=l. 



The momentum equation implied from the impulse approach that 
dH:-* 6 hence Q c = & c 0\Y\h . The definition of the 

coefficient of restitution is independent of position. It may he 
noted that the equations for energy and momentum depend only on the 
inertia distribution, gear ratio and angular velocities. 

In manipulation of these equations it is seen that load inertia 
may be transferred to the£| , shaft by multiplication ^ , the same 

transfer may be accomplished with load friction. Since ^ , and , 
are related to and it may be reasoned that the results 
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for = 1 may be extrapolated to physical systems where /= 1. 

Since a total energy accounting may be made for the system with 
the equations used, the rotational energy lost is attributed to the 
heat of deformation of the gear teeth. This transfer of energy could 
in fact be determined from the equations used. Thus the conservation 
of rotational energy e=l is a special case of a broad interpretation 
of the law of conservation of energy. 
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3. Computer Program Development 
The physical equations of the net systems were programmed for 
solution using the Control Data Corporation 1604 high speed digital 
computer utilizing paper tape program input and magnetic tape output. 
An IBM 717 line printer was used to extract data from the magnetic 
tape. The Control Data Corporation machine library and the U. S. 

Naval Postgraduate School computer subroutine library were used for 
assembly, Runge-Kutta-Gill numerical integration, and decimal output. 

Several changes were made in the forms of the physical equations 
of the net system in order to eliminate duplicate computing operations 
and to fit the equations to a form suited to the variable parameters. 

In the table of symbols, computer mnemonic (m ) terms which are used 
in the assembly subroutine have been indicated by parenthises and will 
be defined when encountered. 

The equation of the system with motor and load combined 

= UV - £u 

was put in the form 

‘ * 7 , • 

- COy\ ■Ac. -_A 6c. wkerc 

A = and 

THETAL = © r ONE = 

THETALD = U - ZETA =? 

* » 

UDOT = 0^ QMEGANSQ 

were the computer m terms used. Underlining will be used in this 
section to denote the transition to the m terms of the computer 
equations. 

The differential equations were solved using the Bunge-Kutta-Gill, 
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(RUNGE), numerical integration method- The RUTJGE subroutine re- 
quired the definition of synonymous terms for the four iterative 
cycles used to produce one extrapolated set of variables. The incre- 
ment of the independent time variable chosen was 0.01 sec. for all 
0.1 and 0.004 sec. for all ^ ^ 6 • | • 

The load free equation 



m # r 

-1- £l B c “ o was put in the form 

* * Lm * 

Be 4- G ° for the computer. 

A separate set of RUNGE synonyms was used for the load free equation 
in order not to destroy the previous computations before they were 
determined of no further value in computing and also in the making 
of program decisions. 



(~)^was designated V and Q was designated VDOT for the computer 



program . 

The equation for the motor alone 

Om + ^ w\ + &c- - JS. 

T ^ 0 'Ad 1 

was put in the form 

^ (&R- &c.) 

d rYl 

- - b. 6 + C ( 6 £ - 



and for RUNGE , the terms 
©fh = THETAM 






YA 



= THETAIM = W 



and 
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YY\ 



= WDOT 



were used. 



The equation representing the law of conservation of momentum and 
the definition of coefficient of restitution were combined to the forms 
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7 XtX — \ P & KVa ( * ^ ^ + -t) 
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Vw 






for the computer program. No provision for additional m terms was 
made to denote the primed values. The additional terms: 



were used. 

Two equations were used to define the boundaries of operation of 
the load and motor when they were acting separately, 



Two major decisions of the computing cycle were the determination of 
the point where the motor and load would float free, FLQATEST , and 
the response of the system to impact of the load and motor gears when 
the boundary conditions on position were met, COMBTEST . 

The first decision, the point of float free, is made on the 
equality or inequality of slopes in the phase plane, i.e., at separation 
N l = Ng. Equating the slopes: 



e = Rest itut p = RH 6 and p z = R Hos$ 



— p yy\ (XV\ d 0 p &YV\ ^ A 




t B =■ 6 



The above equation is solved after computing each point in the 



combined phase trajectory. If 




() the system re- 



mains combined. If 




having 



the response motor alone and load alone immediately thereafter. 



The second major decision, that of whether after impact the 
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response should he that of the combined system or the response of the 
load and motor acting separately, is made on the basis of the resultant 
velocities after impact and upon which side of the backlash the motor 
and load positions are found. 

The computer program, its flow diagram, a description of the pro- 
gram operation, and an explanation of the readout of information with 
examples are presented in the Appendix. 
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4. Methods and Scope of Examination 

Two modes of operation of the computer program were used. The 
first mode considered was the phase trajectory. This mode was used 
to check the operation of the computing cycle. Time, velocity and 
position of the load were printed when the system was comhined. When 
the system was separated, the values of motor velocity and motor 
position were also printed. Although the time increment used for 
computations was either 0.01 sec. or 0.004 sec., printed outputs 
for the first mode were taken only every 0.1 sec. Additional printed 
outputs were taken at the time of contact of the gears just prior to 
and immediately following the solution of the momentum and restitution 
equations. A typical printout of the phase trajectory mode of operation 
is shown in Appendix C. 

The second mode of operation was used for printing only the maxi- 
mum computed positive load overshoot position for each cycle, the 
associated load velocity and the exact problem time of the computation. 

A sufficient number of print-outs for each problem was obtained to 
prove limit cycle existence and average size. Since this investigation 
was concerned primarily with the steady state response, the second mode 
of operation was the one utilized to obtain the majority of the data. 
This read-out method markedly decreased the data reduction time for the 
problem analysis, as opposed to the analysis of steady state conditions 
provided by a full phase trajectory print-out. Since the time required 
for computer read-out was several magnitudes greater than the computing 
time, valuable computer time was saved by this mode of operation. An 
example of the printed output obtained when only maximum load position 
overshoot was of interest is given in Appendix 0 . Approximately 1400 
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phase trajectories were solved using this mode of operation. Each 
solution required an average of three minutes of computer time. 

For the purpose of general examination of the problem it was 
assumed that P = 1, = 1, 0 \\ = f -^L = 1- Specific solutions 

were made with p ^ 1 ( 1 to determine scaling effects. The re- 

sults of these solutions are presented in Section 6, Application of 
Results. 

The program can also be made to solve linear systems by setting 

A = 0. 



The major variable parameters used in this investigation were: 

§ = o.i, 0.2, 0.3, o. 4 , 0.5, 0.6, o.8, l.o 

o , 0.2, 0.4, 0.6, 0.8, l.o 

"^T 

9_5-= 0.1 , 0.2, 0.5 , Q.8 , 0.9 

0 . 9 , 0.8, 0 . 5 , 0.2, 0.1 

A =0.5 radians 



e = 0, 0.6, 0.8, 1.0 
Additional parameters of 



0.04, 0.1, 0.9; and A = 0.01, 0.05, 0.01, 0.15 radians 

Fr 

were used in certain instances to examine particular characteristics of 
the response. 

The program was limited somewhat in that the value of P ^v- = qq 
and zero were exluded due to generation of undefined mathematical 
quantities. 

The value of backlash was made abnormally large, A =0.5 radians, 
to allow an easier interpretation of the non-linear response. Since 
the influence of A was linear, this caused no inaccuracies. All 
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graphs for limit cycle size were plotted with /\ = 0.3 radians. 

The effect of values for /\ other than 0.3 is also shown in graphi- 
cal form. 
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5- Results and Discussion 



Two types of system response were observed. The first type was 
a convergence to zero error, which was a case of static and dynamic 
equilibrium of the system. This is a characteristic of a linear 
second order servomechanism. The second type of system response 
observed was the divergence or convergence to a state of dynamic equi- 
librium. This second type of response is termed a limit cycle and is 
characterized by cyclic travel of the load through the same points 
of the phase trajectory. When the second mode of computer operation 
was used, both the maximum load position per cycle and the period 
of the recorded position were analyzed to determine if the character- 
istic response was a limit cycle. Exact repetition of cyclic values 
was impossible due to the numerical methods used and the computer 
round-off error for the output routine. Typical examples of the 
limit cycle and no limit cycle response are given in Appendices 
D 1 and D 2. 

No case of dynamic instability or divergence of the system with- 
out bound was observed in the investigation. Intuition might lead 
to this same conclusion when the following unique characteristics of 
this system are pointed out: (a) As a maximum limit, energy was 

conserved at the instant of impact, (b) Energy was supplied to the 
system only in finite amounts and at a finite rate by the error 
detector. (c) Energy sinks were present in the system frictional 
elements while energy storage units were present in both the motor 
and load. (d) No undefined limits appeared in the applied physical 
equations, 




O vV\ 



0 and excluded. 
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The results of the investigation are offered in the form of 

charts for parameter areas of limit cycle existence and nonexistence 

in Fig. 1 for e = 0 ., 0 . 6 , 0.8 and Fig. 2 for e = 1 . 0 . The para- 
Z r 2 . 

meters P ti_ and p 3 l are designated abscissa and ordinate 

respectively. The parameter points of examination for which solutions 
were obtained are circled. The zone enclosed by a particular ^ line 
is largely an interpretation by the authors of the results of the in- 
vestigation, and shows the area in which no limit cycle existed. Due 
to the large change of the variable parameter values between solutions, 



the charts are not exact. All circled intersections interior to the 



s 



line indicate parameters which resulted in no limit cycles, e.g.. 
Fig. 1 for e = 0, 0.6, 0.8 at a ^ = 0.4, the values of = 



Ft 



0.6, 0.8 at a ^ =9*0 and the value of ^ it = 0.6 at a 

0 wt 

= 4.0 resulted in no limit cycle. Similarly, for = 0.5, 



3 Wt 



the values of = 0.4, 0.6, and a O =1.0 as well as all 

TV 

other parameter intersections interior to the =0.5 curves, re- 

sulted in no limit cycles. 

All circled parameter intersections exterior to particular ^ lines 
define paramters which were examined and resulted in a steady state 
limit cycle. The area to the right of the single ^ =1.0 line 
designates the parameters for which no limit cycles were observed. 

Limit cycles were observed for all ^ ^ 0.5 for the parameters examined. 

The investigation was made with ^ =0.5 radians; however, re- 

sults of the investigation of variation of limit cycle size indicate 
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that existence of the limit cycle is independent of the size of back- 
lash in a nonlinear system for a A. greater than zero. Variation of 
A had the effect of a linear variation of limit cycle size only. 

The zones depicted on the existence charts are approximately 
represented by straight lines; however there is reason to believe that 
a complete examination would show a curvature in the lines. This 
observation will be elaborated upon later when the relationship between 
the existence charts and the figures of limit cycle size is discussed. 

From Fig. 1 it is seen that there is only a very small difference 



in the limit cycle existence zones for plastic impact (e = 0) and the 
intermediate values of the coefficient of restitution, e = 0.6 and 0.8. 
When comparing the elastic impact (e = l) in Fig. 2 to the plastic and 
intermediate cases in Fig. 1, the most obvious point of comparison is 
the drastic reduction of the areas in which limit cycles did not result 
for e = 1 . 0 . However, it may be generally observed that with any e 
held constant, the area in which a limit cycle could not be obtained 
increased with the increase of system ^ . 

For all values of e, there is an apparent common enclosed zone 

centroid which is described by the system parameters yy f ^yF^ =0.5 

Vi- 
and p ^L/^m = 1*0' These parameters result in the system time con- 
stants = si nce £0 n = 1 a nd = 1 were chosen for the 

investigation. There is also an apparent symmetry about the centroid 
of opposite quadrants which can be related to the fact that since 
60^= 1 and p = 1, the points diametrically opposite across the 
centroid of the figure (using the decimal values of the parameter inter- 
section) represent an exact exchange of component time constants. (See 
Fig. 2 for example). Two diametrically opposite points may have 
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exactly the same motor and load time constants; however, they are 
descriptive of two different systems since the friction and inertia 
distributions are different for the two points. 

Figs. 5 through 26 present the maximum positive error in radians 
of the limit cycle for a unit step input and with backlash equal to 
0.3 radians. The use of the charts is explained in Section 6, Appli- 

tL. 



cation of Results. Each figure is presented with p ^L/^T an< ^ 6 



constant. The independent parameter 



,J is indicated as abscissa 



? V- 

and curves of constant 3 are plotted with limit cycle error as the 



ordinate. The results were presented on a log log plot since the 
observed limit cycle errors encompassed a range from 0.01 to 1.0 
radians, and the appearance and disappearance of limit cycles and 
the change of limit cycle error was more readily apparent on the log 
log form. 

In general as e increases, the size of limit cycle decreases, all 

other parameters remaining constant. However this is not without a 

few exceptions; the generalization is most accurate for large size limit 

cycles. It may be noted that when viewing the individual charts, if a 

k 

limit cycle exists at very low P J /j ratio there is a tendency for 

\ L M 



the same system to have a smaller limit cycle for higher 



ratio . 



J/J 

\ L m 

In addition, for these same conditions, if limit cycles exist for 

several values of ^ , there is a convergence toward the same size of 

limit cycles for each ^ as pj/j increases. At the value 

\ L m 

P 



f /F = 0, a limit cycle will always occur for any value of e or 
L T 



s- 

Viewing the charts in sequence, with either e constant or <p f^/F 
constant, it is noted that there is continuity of pattern flow between 
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adjacent charts. This continuity between charts appears to form a 

symmetrical pattern of constant lines which can be interpolated for 

L 

intermediate values. When the figures for p f /F = 0.4 and 0.6 

V L' t 

are analyzed, this symmetrical pattern is observed in each chart about 

p J /J a l. The limit cycle decreases and disappears at the 
V L m t 

The charts of f)f/F = 0.4 and 0.6 
V 1/ T 



opposite extremes of p J /j 

1 L m 

are symmetrical to each other, while f^/F^ = 0*2 and 0.8 are 
similar. At the extreme values of ^ f /F^, = 0 and 1.0, no symmetrical 
comparison can be made. An additional point of symmetry about 
^ f^/F^ ~ can be visualized for all values of e. These para- 

meters for symmetry are noted in the discussion of the existence charts . 
The irregularities in the curves of medium values of ^ can not be 
explained by the authors; however, it may be mentioned again that the 
points of symmetry are characterized by exact exchange of component 
time constants. 

The noted symmetry of the limit cycle size charts shows a direct 
relationship to the symmetry found in the existence charts. Since the 
existence charts show contours of constant ^ for a surface of zero 
limit cycle size, one is led to the possibility of a three dimensional 
presentation to show a variation of limit cycle size for the same para- 
meters used in the existence charts. Such a three dimensional figure 
for constant e using ^ f^/F^, and ^ as the axes is shown 

in Fig. 27 ♦ Other surfaces could be added above this zero limit cycle 
"floor" to show the variation of a given limit cycle size by contours 
of ^ . The ultimate result of this line of investigation would be 
the evaluation of the poles and zeros for a confoimal mapping presen- 
tation. Unfortunately, time did not permit the authors to investigate 
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this avenue flirt her. 

The results of the investigation of the dependence of limit cycle 
size on the amount of backlash are presented in Figs. 28 and 29 using 
coefficients of restitution of 0 and 1.0 at a ^ =0.6. A thorough 
investigation was made at this ^ since it is an average system 
parameter often used in design work. The limit cycle size was found 
to vary directly with backlash size without exception between the 
values of 0.3 and 0.01 radians. When limit qycles existed at ^ = 0.3, 
there was no disappearance of limit cycle with decreasing magnitude of 
A. . Where no limit cycle existed, at /\ = 0.3, limit cycles con- 
tinued to be non-existent throughout the range of examined. 

Several solutions of problems with = 0.6 and e = 0.6 &0.8 and = 
0.1&e=1.0&0 with various values, proved the linear variation 
of limit cycle size with backlash size for all values of e. 

When a servomechanism under steady state conditions is disturbed 
by a small perturbation, it will return to the original steady state 
conditions. However, the transient response due to the perturbation 
may be quite different from the transient response which occured due to 
the original signal input. The concepts of energy and static and dy- 
namic equilibrium mentioned in Section 1 may be better understood as a 
result of the perturbation technique if the usual initial conditions of 
A (0) = A (0) = 0 are considered as a disturbance from the stability 
conditions predetermined by selection of system parameters and signal 
input. 

This same small perturbation technique was simulated mathematically 
using the digital computer to obtain several solutions. The cases 
examined and their initial conditions were: 
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a) ^ = 0.6, ^J L /J m = 0.8/0.2, p"f L /F T = 0, e = 0, A=0.3 

-6 c c (0) = 0 ; -6 c (0) = 1.2 

4 m (0) = 0 ; 0m (0) = 1.0 



system separated. 




system combined 



c) '§=0.4, P'j/J = 0.5/0. 5, 

\ L m 

c (0) = 0 j & c (0) = 1.005 

* 

■0- m (0) = 0 ^-m = 1.00 

•J 




system separated 

In the three cases examined, the same limit cycle was obtained, 
both in size and period, as that resulting from the normally used 



between the perturbation response and the original limit cycle size 



response were observed. These were convergence to limit cycle in the 
first two cases and divergence to limit cycle in the last case. 

The results obtained by the small perturbation technique supported 
the validity of the analysis of the problem by digital methods. This 
result also concurs with the findings for the case of e = 0 by Knoll 
and Narud, Ref. d, in that the steady state limit cycle size is completely 
specified by the choice of system parameters and independent of system 
initial conditions. 

Mention is made here that the computer program could not be run 
normally for the perturbation method. The computer was stopped after 




for the three cases was 18 x 10 - ^ radians. Two types of transient 
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the print out of system constants, then the initial conditions of the 
servo system were inserted manually. The program was then restarted 

at the desired point for computing, either as a combined system or 

with the motor and load separated. The program was then allowed to 

operate to the normal completion of the solution. 

The work of Ref. d reported that the system transient response 
for the case e = 0 arrived at a steady state limit cycle by two 
possible phase trajectories. One of these phase trajectories was 
smooth convergence to the limit cycle, while the second was a converging 
overshoot toward the inside of the limit cycle followed by a divergence 
to the limit cycle. From the results of this investigation for the 
cases e = 0, 0.6, 0.8 and 1.0, similar transient indications were ob- 
tained. 

Several computer solutions were accomplished for the cases =0.3 
radians, e = 0, 1.0, D = 1.0, It) n = 1.0, K = J / 1.0. In addition 



purpose of these solutions were: 

1. To determine the general applicability of the computer program. 

2. To verify the scaling developed in Section 2, Development of Equations 

3. To determine the applicability of the results of this thesis. 

The following results were determined: 

A. Two solutions were obtained (e = 0 & 1.0) for the case ^ = 0A, 



The results of these solutions were the same (time, load velocity, load 



various combinations 




obtained. The 




position, and limit cycle size) as the solutions obtained for 





nondimens ional parameters (j) n = 1.0, 






^ = 1.0, J = 1.0, K = 1.0. 
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B. Solutions were obtained for e = 0 & 1.0 for the case = 0.4, 

= 1.0, 00 n = 2.0, J = 1.0, K = 4.0, = 0.4, J?Ti_ = 4.0. 

Ft ' Twi 

The results of these solutions yielded the same size maximum overshoots 



e 



§ = 0.4, 



and limit cycle size as the nondimensional solutions for 

Pjh_= 0.4, ^ .0 1 - 4.0. The times of maximum overshoots in 

' Ft 3" vv\ 

this case were one~half the times of occurence of the same overshoots 



for the nondimensional solution. The velocity at a maximum overshoot 



was twice the value of velocity obtained for the same overshoot in the 
nondimensional solution. This same comparison can be made for Cases A 
and B. The results of these bnd other solutions prove the validity 



formulae used in Section 2, Development of Equations, concerning the 
transformation to the coordinate system. (*) 

C. Solutions were obtained for the case e = 0 & 1.0, 'S = 0.4, p = 0.5, 

£l)n = 1.0, 4j_ = 0.8, Tl_ = 4.0, = 0.2, =1.0. The 

f~T T Ft t wv 

results of these solutions were the same (time, load velocity, maximum 



overshoots and limit cycle size) as for the nondimensional solutions 

= o.4, P = 1.0, 6d n = 1.0, = 0.2, P 1 _7u= 1.0. 

v Ft ' 

D. Solutions were obtained for the case e = 0 & 1.0, ^ = 0.4, 

P = 0.5, lj) n = 2.0, = 0.8, = 4.0, ( £lt- =0.2, ^ Jl = 1.0. 

Pf TTyv'x FT "3 yv\ 

The results of these solutions yielded the same size maximum overshoots 



and limit cycle size as the nondimensional solutions for = 0.4, 



& 



) ^L= 0.2, P = 1.0. The times and velocities in this case were 

Ft , . Sw 
scaled by 60 n as in Case B. 



Several phase trajectories using e = 0 were compared with those 



obtained by New, Ref. e. Generally good correlation was obtained on 
the value of limit cycle error and in the transient response. Small 



disagreements were expected, resulting primarily from differences in the 
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programming procedures and analysis of results by read-out methods. 



Sample phase trajectories for the cases of e = 0 & 1.0 are pre- 
sented in Figs. 30 and 31 . Although the transient analysis of the 
system is beyond the scope of this work, these figures may aid in 
understanding the physical problem. 

Of academic interest was the accidental operation of the computer 



observed. One was a nearly pure oscillatory system from the first 
overshoot and thereafter another was a very slow convergence or slow 
divergence of the oscillations after the first overshoot. The con- 
vergence or divergence probably resulted from errors generated in 
numerical iterations and approximations. 



program for several solutions using 




Several responses were 
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6. Application of Results 

In this section, the results of this thesis will be applied to 
hypothetical physical systems. Methods will be indicated by which physical 
systems may be analyzed or designed. Attention is called to the assump= 
tions stated in Section 2 under which this study was made. 

It is desirable to find an equation relating certain system 
variables, therefore the equations: 



a) 



U)^ = P ^ 

vi yyy -V T L 

a v\ d 



- P- 



k 






( 2 ) 



t j3 f 

P ; 






( m v n 
P^L J 



P 









will be manipulated to the form 

(5) ^ = -1 • 

0 w\ 
or 



+ P 3l 



DJl 

T*1 



M 






■h \ 



4- ^ 

1 ml 



VTtn \| \ + p"2t 
vHT \nr ' dm 





C - J- £ , \ f "- 1 4-1 


1 -1 

P 


00 




- 3 vn C 4- ^ vW ^ 






In view of the form of Equation (4) for ^ and recalling that the 

variable parameters used on the figures were P 3-b and P ^i=. = 

^ r \ ' 3 "™ 

V •“ L. = r .. , \ it is seen that the expression for -P 

r *■ 1 O 

T \ P V U p 1 - r 1- 

in fact includes the other two variable parameters. 
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Fran discussion concerning the effects system natural frequency, 
it will be recalled that the value of overshoot and hence limit cycle 
are not affected by CO n; however, the time and velocity are scaled 
for n ^ 1.0. Thus the transient response characteristics may 

be related to the same parameters used in these applications. 

If a motor and load are selected for a system and the limitation 
is imposed that system operation shall not result in a limit cycle, 
Equation (4) and the limit cycle existence charts may be used to 
determine the unspecified system parameters , K and ^ . The 

coefficient of restitution for the proposed gear train material may 
be determined by the use of the equation in Section 2. 

To proceed, select a value of K and (p , then determine the value 
of ^ from Equation (4). Enter the limit cycle existence chart for 

the value of e, system and the parameters P ^ *- and — — ^ . 

-*-1 



p z (i 



Determine if the selected parameters describe a point of no limit 
cycle. Since limit cycle existence is independent of backlash ( f 0 
excluded), backlash is not a variable parameter. Successive trial 
values of the unspecified parameters may be required to meet the 
limitation of no limit cycle. If it is also required that certain 
transient response characteristics be met, these conditions could be 
examined by a similar method. 

In the event that no acceptable parameters are found to satisfy 
the system requirements for no limit cycle, the figures for limit cycle 
size may be examined to obtain a minimum size limit cycle. The method 
of solution is similar to that previously described. 

Example: The material intended for use in the gears has an e = 0.6. 
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The motor and load have the parameters fin = 0.64, Jm = 0.25, f T = 0.64, 

.L 

= 1.0 and it desired that the system operate with no limit cycle. 

If the values are selected K = 4.0, ^ = 0-5, the values = 1.0 

4m =4.0. L ~ ~ -■ ^ ™ 



and i M 



e' f 



= 0.2 result. 



. 5 






f ^ +■ 1 

Solve Equation (4) for t - .5(.64)(5+l) 

= 0.4. Enter the chart for limit cycle existence. Fig. 1, 
e » 0.6, 'C - 0.4, — - J - 0.2 . C? 31 • 1.0 

* Sin-' ' ~ 



, thus 

2. 



p*C< 






It can be seen that limi-t cycle will exist for these parameters. 
Since K is most easily varied in equation (4), a new K is selected 
at 1.0 in an effort to find parameters for no limit cycle. It is seen 



that 



t 



K = U 



t 



K-- H.o 



= 0,4 (?) 



with all other parameters constant, thus ^ =0.8 Enter Fig. 1 

rtth e=0.6 <).* ^ , J () J 



0 wi 



J 



- 0,1 



fi 



_J21 + l 

These parameters describe system operation with no limit cycle and are 
well within the area limits. The transient response of this system 
may be undesirable and various combinations of p and K would have to 
be tested to satisfy the additional requirements of the problem. It 
will be restated that the existence areas of Figs. 1 and 2 are not exact. 

The figures of limit cycle size may be used in the same manner as 
the existence charts. The same parameters of the first case are 
assumed: K = 4.0, p = 0.5> ^ = 1.0, 

O VY\ 
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ru 



= 0.2 



K - 0.4 



e = 0.6 A =0.1 radians. Enter Fig. 10 with the parameters given. 

The ordinate, Limit Cycle Error ( 0. 3 \ = 0.13 radians, is obtained. 

V A / 

The maximum positive error in the limit cycle = 0.13 (O.l) = 0.0^3^ 

radians. Enter the same figure for e = 0.6, ^ =0.8 ) ( Ju = 1.0 

l „ “ 

i » = 0.2 , = 0.5 radians; this is the second case of exami- 

pFu 1 ; 

nation for existence. No \T = 0.8 curve exists; hence no limit cycle 
exists for this choice of parameters for any nominal value of backlash. 



) 
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7. Conclusions 



Within the limitations of the assumptions made and the scope of 
this thesis, it is concluded that: 

A. Two types of system steady state response result. The type of re- 
sponse is independent of backlash size, system natural frequency and 
initial conditions. The types of response are convergence to signaled 
input position and convergence to a limit cycle. 

B. If a limit cycle exists in the steady state response of a system, 
the size of the limit cycle varies directly with the size of backlash, 
all other parameters constant*. 

C. The greatest change in system steady state response occurred for 
the case when system rotational energy is conserved (e = 1.0). 

D. A point of symmetry for the type of system steady state response 
exists with equal motor and load time constants. For medium to high 
values of system damping coefficient, symmetry with respect to limit 
cycle size exists for an exchange of motor and load time constants. 

E. The results of this thesis are presented in a form which may aid 
in system design or analysis. 
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8. Recommendations 



From the results of this thesis it can be seen that there remain 
many fertile areas for further investigation. For studies of this 
type it is strongly recommended that the digital computer be used 
because of its speed and versatility. The limits of investigation 
depend to a great extent on the ingenuity of the programmer. Some 
recommended, areas of future investigation are: 

a) More precise definition of the limit cycle existence zones 

and a more thorough coverage of the limit cycle size of similar systems 
having the same type of parameters. A three-dimensional presentation 
as mentioned in the discussion could be obtained. 

b) Include torsion and other deformations of the mechanical system 
in the physical equations of the system to determine their affect on 
the system response. 

c) An analysis of the limit cycle existence zone and transient 
response of a system to determine if there is a relationship between 
the two. An analytical expression for the possibility of limit cycle 
existence might be developed in terms of system parameters. 

d) Examine the transient and frequency response of the system due 



to the computation section of the existing computer program. 

e) Since the physically realizable quantities of position, velocity 
and acceleration of both the motor and load are available as computed 
quantities, the study of the compensation problem and its optimization 
are obvious areas for further consideration. 

f) The general study of other higher order nonlinear systems is 
suitable for digital computer programming since most numerical inte- 




input functions. This would require minor modifications 
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gration methods will handle systems having equations of any order. 
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APPENDIX A 

COMP U TA T/OSV FLOW £>/A G-PA M 

for 

CDC /604 D/G-/TAL COMPUTER 




7Z 











from MOTOR. 3Q 




73 




from S9, GO 




from S£,6 





from 4*3 




83 



EXIT’ _ 

84 





ref urn to progrmtr 
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APPENDIX A (continued) 



1. 



2 . 



3- 

k. 

5- 

6 . 



7- 

8 . 

9- 



10 . 



ll. 



12 . 



13. 

14 . 



EXPLANATION OF FLOW DIAGRAM BLOCK SYMBOLS 
Clear fault stop. 

Print last overshoot stored in LIMIT3UF (symbol 888 ). 

Print out position in phase trajectory of fault stop 

(symbol 666 ). 

End of file. 

Actuate fault stop. 

Load zeta from Index 1 (b ). 

1 

Print last overshoot stored on LIMITBUF (symbol 888 ) . 

Print position in phase trajectory where it was manually 
stopped (symbol 1961). 

Stop if Zeta - 0. 

2JO). - a 

Is zeta - 0,1? 

Is zeta less than 0.1? 

Set up computation time of 0.01 seconds for problem and 
set up for print at every 0.1 seconds when zeta is greater 
than 0.1. 

Set up computation time of 0.00^ seconds for problem and. set 
up for print at every 0.1 seconds when zeta is equal to or 
less than 0.1. 

Set up limit cycle print-out routine for 20 overshoots and 
average the last 8. 

Set up number of limit cycle print-outs and averaging routine 
depending upon the size of zeta: 



75 



f 


Overshoots 


Start 

Averaging 


Divide 
by * 


Comp. 

Time 


^ 0.8 


6 


4 


2.0 


0.01 


0.6 


7 


5 


2.0 


0.01 


0.5 


8 


5 


3.0 


0.01 


0.4 


10 


7 


3-0 


0.01 


0.3 


14 


10 


4.0 


0.01 


0.2 


20 


12 


O 

CO 


0.01 


< 0.1 


17 


12 


5.0 


0.004 



*The last row printed has this number printed on the far left 
if the problem has a normal solution. This number indicates 
how many overshoots were averaged. The average value has 1.0 
subtracted from it and is printed in the far right of the same 
last row. 

15 . Page re j ect . 

1 6 . Clear buffers and set up for new problem. 

17 * Restitution, e, from Index 3 (b^,). 

Backlash, , from Index 4 (b^). 

J rr, from Index 5 (b,_ ) • 

5 

Jt “ Jw K - Jtn 

jjt- from Index 6 (bg). 

(y fT y - f*- j f !yf L - & 
v° 2fL = & j - 0 

A r 

18 . Have all values been cycled through? 

19. Set STOPNOW after this run to stop absolutely. 



= 0 ? 



20. Backlash, ^ , 



21. Print consecutive number of this parameter run and all the 
parameters used: 

£j°6J»=A J'Yjl Vf t & 

J L JmV 9 *! frn f L /° 2 

k K /j^c Vj.-o / 

22. Set program so always acts as combined system. 

23. Print consecutive number of this parameter run and selected 
parameters for linear system: 

Zf6>*=A K GO*? 



/° /° 



2k. Solve the second order equation by Runge-Kutta-Gill numerical 

integration every 0.01 seconds problem time when zeta is greater 

than 0.1, otherwise every 0.004 seconds. Initial conditions 

• •« 

are 1«0, © c = ~ ® an< ^ = ° r P re “ 

viously computed point. Next point is computed by four 
iterations of: 

0 C ~ ^ h (©r ~ ©c) “ & jf* CO h 0c 

25. Set EXIT from print routine. 

26. Load time, a c and in print buffer. 

27. Go to LDdTCYC to print only the overshoot. Otherwise print 
the phase trajectory points for each 0.1 seconds up to a maxi- 
mum of 12 cycles. 

28. Store 0 for and © m 
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29. 



30. 



31. 



32. 



33. 



34 . 



35- 

36. 



37- 



38. 



39. 

40. 

41. 

42. 

43. 

44. 
45- 
46. 



Is problem time 0.1 sec.? 

Is problem time greater than 0.1 sec.? 

o 

Print the number of point computed, time, © c and ©c . 
EXIT (automatically set for desired jump-out.) 

O® ° 



Store time, ©C ’ e c and © c for motor and load initial 

conditions. 

Is slope N t = N_? 

-L' S 

(©R-G>c)d0n_ 3 __ * 

sr- 1 — eLy ^>v\T- J-h- = positive value? 

J L 

Is Q negative? 

Store © c - A irr 



Solve 



©c = -fycK&c) 



for ©c and e c 



by Runge-Kutta-Gill. 




for and 0* by Runge-Kutta-Gill 

0 L ” ^ o ? 

0 L ~ /° Q A = negative? 



Set EXIT from print routine. 

Set EXIT from print routine. 

Load time, and Q c , 

o 

Load and - 

Print number of point computed and stored values. 
EXIT (automatically set for desired jump-out). 
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4-7 • Solve 

© / 

Oc = 



I KKV 






<9 m (i + e) -t © c (/? -A - e) 



©m 



©* s 



& - e(/<?erh +■ ©c) 




o / 

-h 



7- (©* - ©«) 



4-8. Set EXIT from print routine. 

© 

4-9. Is ©c positive? 

© 

50. Set up to count one phase trajectory when ©c is negative. 

51 . is 0 C _> l.o? 

52. Have 12 cycles of phase trajectory "been completed? 

53* Arrange for no counting until &c is positive again. 

* A 

5^. Is ©c negative? 

55 • Set up — (jr ® “ S>c) - 

* © 

56 . Is "* greater than preselected £ ? 

57 * Set up combined system. 

58 . Is © c positive? 

© 

59. Is ©c negative? 

60. Is G> c - l positive? 

6 1. Is — I negative? 

62. Is 0c“ I negative? 

6 3- © c -^0m =o? 

64. O c -fBn=0l 

65 . Is <3^ positive? 

66. Is $ c ~yO Prepositive? 
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6 7- 



68 . 

69- 

70. 

71. 

72 . 



73. 

74 . 



75- 

76. 



77- 



78. 



79- 

80. 

81. 

82. 

83. 

84. 

85. 

86 . 
87. 



Is Slope = Kg 

C I - Sc) U)n 



O 

©C 



- + 



ft 



Modify FLOATEST to keep system combined. 



when ©e is negative. 

Is 0^ — positive? 

Store ©L. ~ A in @ ^ 

Store Oc in 0 ^ 



- positive value? 

Takes off modification 



Is ©C ^ I negative? 

0c = 0 ? 

0c " 0? 

Is ©c ~ ©^positive? 

Is @ c -y/’S ©** positive? 

Modify FLOATEST to keep system combined. Takes off modification 

c 

when is positive. 

Have ll/2 minutes of real time elapsed since the last print-out? 
Is 0 C < 1.0? 

Set up to print when in the fourth quadrant. 

Check if in the fourth quadrant. 

Is the new ©e greater than the previous 6>c ? 

* 

Store the new time, &c and ©c 

EXIT (automatically set for desired jump-out). 

Print the last overshoot stored in LIMITBUF (symbol 888 ) 

Stop real time clock. 

Print position in phase trajectory where it was stopped by the 
real time clock (symbol 999 ) • 
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88. 

89. 

90. 

91. 

92. 

93. 

94 . 
95- 
96. 
97- 
98. 

99- 

100. 

101 . 
102 . 

103. 

104 . 



Return to "beginning to start new problem if new parameters 
remain to be evaluated. Unconditional stop if all parameters 
have been evaluated. 

Has the required number of overshoots been evaluated prior to 
starting the averaging routine? 

Add the just evaluated to the previous sum for averaging. 
Add one to the count prior to averaging. 

Clear time clock. 

o 

Print time, and a c . 

Required number of total print-outs? 

Clear buffers. 

Obtain the average of ©c- 
Subtract 1.0. 

Print the average of 1.0 and how many overshoots were 

utilized to obtain the average. 

Is ©c less than 1.005 radians? 

Clear buffers. 

Start real time clock. 

EXIT (automatically set to desired jump-out). 

Has 30 seconds of real time elapsed since the last print-out? 
Print the position in the phase trajectory where it was stopped 
by the real time clock and being less than 1.005 radians 

(symbol 1005) . 
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Jump Switches: 



- For LIMITCYC , allows the maximum overshoot print- out only. 
With down, a phase trajectory is printed for a maximum 
of 12 cycles. 

- Set the number of print-outs depending on zeta. Change 
the number of overshoot values used to obtain the average 
value. 

J - Use the real time clock for an automatic recycle to a 

3 

new parameter when have: 

a. ll/2 minutes maximum after the last print-out. 

b. 30 seconds maximum and is 1*005 radians or less. 



Stop Switches: 

S^ - Stops at the end of this parameter run or stops at the end 
of the manual print-out. 
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APPENDIX B 



COMPUTER PROGRAM 









00070 


ORG 




00007 


00007 


74 


0 


FAULT EXF 


0 


00070 




75 


0 


40701 


SLJ 


0 


FAULPRIN 










ORG 




40000 


40000 


74 


0 


00100 


START EXF 


0 


00100 


40001 


75 


0 


40004 


SLJ 


0 


MEW Z ETA- 1 


12 


1 


40623 


LA2Y2ETA LDA 


1 


ZSTAINDX 




20 


0 


40425 


sta 


0 


ZETA 


40002 


50 


1 


00000 


ENI 


1 


0 


40003 


75 


0 


40000 


SLJ 


• 0 


START 


75 


0 


40674 


MANSTART SLJ 


0 


MANUAL 




50' 


0 


00000 


EMI 


0 


0 


4 OOO 4 


12 


0 


40565 


LDA 


0 


NEWSTOP 




20 


0 


40563 


STA 


0 


STOPNOW 


40005 


12 


0 


40425 


NEWZETA LDA 


0 


ZETA 




22 


0 


40564 


AJP 


0 


STOP I NDX 


40006 


32 


0 


40665 


FMU 


0 


TWO 




32 


0 


40426 


FMU 


0 


OMEGAN 


40007 


20 


0 


40404 


sta 


0 


A 




12 


0 


40624 


LDA 


0 


ZETAINDX+1 


40010 


31 


0 


40425 


FSB 


0 


ZETA 


40011 


22 


0 


40433 


AJP 


0 


OUIKCOMP 


22 


2 


4Q433 


Ao? 


2 


QUIKCOMP 




12 


0 


40402 


LDA 


0 


POINT01 


40012 


20 


0 


40154 


STA 


0 


TABLES-i-l 




20 


0 


40173 


STA 


0 


TABLEL+1 


40013 


20 


0 


40207 


STA 


0 


TABLEM+1 


12 


0 


40661 


LDA 


0 


COUNT 10 


40014 


20 


0 


404 00 


STA 


0 


TENTHSEC 


50 


0 


00000 


ENI 


0 


0 


40015 


75 


2 


40437 


ZZ SLJ 


2 


BIGZETA 




50 


0 


00000 


EMI 


0 


0 


40016 


12 


0 


40652 

40615 


LDA 


0 


STOP 20 




20 


0 


STA 

LDA 


0 


STOPZETA 


40017 


12 


0 


406o2 


0 


COUNT 12 




20 


0 


40617 


sta 


0 


COUNTZET 


40020 


12 


0 


4067? 


LDA 


0 


EIGHT 




20 


0 


40620 


STA 


0 


DIVIDE 


40021 


10 


0 


00010 


ENA 


0 


10 


61 


0 


40562 


SAL 


0 


STOPRINT 


40022 


75 


A 


71000 


GOAROUND SLJ 


lx 


DECOF 


50 


6 


00000 


ENI 


6 


0 


40023 


02 


0 


40427 


02 


0 


FORZEROS 




04 


0 


00000 


04 


0 


0 


40024 


12 


0 


40664 


LDA 


0 


ONE 


20 


0 


40156 


STA 


0 


UDOT 


40025 


10 


0 


00000 


ENA 


0 


0 




20 


0 


00000 


STA 


0 


0 


40026 


20 


0 


40162 


STA 


0 


THETA 


20 


0 


40161 


STA 


0 


THETADOT 


40027 


20 


0 


40157 


STA 


0 


U 


20 


0 


40155 


STA 


0 


1 


40030 


20 


0 


40377 


STA 


0 


INDEX 


20 


0 


40615 


STA 


0 


' NDEXAVE 


40031 


20 


0 


40351 


STA 


0 


PRI NTBUF-r3 




20 


0 


40352 


STA 


0 


PRINTBUF+4 


40032 


20 


0 


40353 


STA 


0 


PRI NTBUF+5 


20 


0 


40614 


STA 


0 


LIMINDEX 


40033 


20 


0 


40607 


STA 


0 


LIMITBUF+2 




20 


0 


40610 


S l A 


0 


LIMITBUF+3 


40034 


20 


0 


40611 


sta 


0 


LI M I TBUF+4 


20 


0 


40612 


STA 


0 


LIMITBUF+5 


40035 


20 


0 


40613- 


sta 


0 


L I M I TBUF+6 




12 


0 


40621 


LDA 


0 


PRINTSET 


40036 


20 


0 


40326 


STA 


0 


0K2PRINT+1 


20 


0 


403 35 


STA 


0 


B0UNPRIN+4 


40037 


12 


0 


40622 


LDA 


0 


NEWLIMPT 




20 


0 


40544 


STA 


0 


LIMPRINT+1 


4o040 


12 


0 


40253 


LDA 


0 


FLOTDAT A . 


- 


20 


0 


40071 


STA 


0 


FLOATEST 


40041 


12 


0 


40254 


LDA 


0 


FLOTDATA-i-1 


20 


0 


40072 


sta 


0 


FLOATEST+1 


40042 


75 


0 


40452 


NEWINPUT SLJ 


0 


I NDXVALU 


50 


0 


00000 


ENT 


0 


00000 



4004-3 


12 


0 


404io 




LDA 


Q 


DELTA 


40044 


22 


0 

0 


40500 




AJ? 


0 


LINEAR 


12 


40513 




LDA 


0 


NOSINK 


40045 


20 

75 


0 

4 


Toll 




S i A 
SLJ 


0 

4 


FLOATEST 

DECOF 




50 


0 


00000 




ENI 


0 


0 


40046 


01 


0 


40404 




01 


0 


A 




06 


0 


00001 




06 


0 


1 


40047 


72 


0 


40046 




RAO 


0 


7-1 


40050 


2? 


§ 


71000 

40412 




S W 


4 

0 


DECOF 
MO INERT 




06 


0 


00000 




06 


0 


0 


40051 


72 


0 


40510 




RAO 


0 


LINEAR+10 




75 


4 


71000 




SLJ 


4 


DECOF 


40052 


01 


0 


40420 




01 


0 


KMOTCONS 




06 


0 


00000 




06 


0 


0 


40053 


10 


0 


00000 


combined 


ENA 


0 


0 


40054 


20 


0 


40160 


STA 


0 


QX 


20 


0 


40163 




STA 


0 


QY 




75 


4 


60200 




SLJ 


4 


RUNGE 


40055 


00 


0 


40153 




. 0 


0 


tables 




00 


0 


40164 




0 


0 


DERIVES 


40056 


75 


4 


60201 




SLJ 


4 


RUNGE+1 




50 


0 


00000 




ENI 


0 


0 


40057 


75 


0 


40310 




SLJ 


0 


COMB PR IN 


50 


0 


00000 




ENI 


0 


0 


40060 


12 


0 


40155 




LDA 


0 


T 




20 


0 


40174 




STA 


0 


TL 


40061 


20 


0 


40210 




STA 


0 


tm 


12 


0 


40156 




LDA 


0 


UDOT 


40062 


20 


0 


40175 




STA 


0 


VDOT 




33 


0 


40411 




FDV 


0 


RHO 


40063 


20 


0 


40211 




STA 


0 


W DO T 


12 


0 


40157 




LDA 


0 


U 


40064 


20 


0 


40176 




STA 


0 


V 


40065 


20 


0 


40200 




STA 


0 


THETADL 


33 


0 


40411 




FDV 


0 


RHO 




20 


0 


40212 




STA 


0 


W 


40066 


20 


0 


40214 




STA 


0 


THETADM - 




12 


0 


40162 




LDA 


0 


THETA 


40067 


20 


0 


40201 




STA 


0 


THETAL 




33 


0 


40411 




FDV 


0 


RHO 


40070 


20 


0 


40215 




STA 


0 


THETAM 




50 


0 


00000 




ENI 


0 


0 


40071 


13 


0 


40162 


FLOATEST 


LAC 


0 


THETA 


30 


0 


40664 




FAD 


0 


ONE 


40072 


32 


0 


40424 




FMU 


0 


OMEGANSO 




33 


0 


40161 




FDV 


0 


THETADOT 


40073 


31 


0 


40404 




FSB 


0 


A 


30 


0 


40423 




FAD 


0 


B 


40074 


22 


2 


40053 




AJP 


2 


COMBINED 




12 


0 


40161 




LDA 


0 


THETADOT 


40075 


22 


3 


40270 




AJP 


3 


MODIFY 


50 


0 


00000 




ENI 


0 


0 


40076 


10 


0 


00000 


LOAD 


ENA 


0 


0 




20 


0 


40177 




STA 


0 


QL 


40077 


20 


0 


40202 




STA 


0 


OLL 


20 


0 


40213 




STA 


0 


OM 


40100 


20 


0 


40216 




STA 


0 


QMM 


75 


4 


60200 




SLJ 


4 


RUNGE 


40101 


00 


0 


40172 




0 


0 


TABLEL 


00 


0 


40203 




0 


0 


DERIVL 


40102 


75 


4 


60201 




SLJ 


£ 


RUNGE+1 


50 


0 


00000 




ENI 


6 


0 


40103 


. 75 


4 


60200 


MOTOR 


SLJ 


4 


RUNGE 


50 


0 


00000 




ENI 


0 


0 


£0104 


00 


0 


40206 




0 


0 


TABLEM 




00 


0 


40217 




0 


0 


DERIVM 


40105 


75 


4 


60201 




SLJ 


4 


RUNGE+1 


50 


0 


00000 




ENI • 


0 


0 


40106 


13 


0 


40215 


CONTEST 


LAC 


0 


THETAM 




32 


0 


40411 




FMU 


0 


RHO 


40107 


30 


0 


40201 




FAD 


0 


THETAL 


22 


0 


40112 




AJP 


0 


BOUNCE 



84 



40110 


22 


3 


40112 




A JP 


3 


BOUNCE 




31 


0 


40410 




FSB 


0 


DELTA 


40111 


22 


3 


40317 




A JP 


3 


FLOTPRIN 




50 


0 


00000 




ENI 


0 


0 


40112 


75 


0 


40331 


BOUNCE 


SLJ 


0 


BOUNPRIN 




50 


0 


00000 




ENI 


0 


0 


40113 


12 


0 


40200 




Lda 


0 


THETADL 




20 


0 


40212 




sta 


0 


W 


40114 


12 


0 


40417 




LDA 


0 


RHOSQ 




33 


0 


40405 




FDV 


0 


INERTRAT 


40115 


31 


0 


40407 




FSB 


0 


RESTITUT- 




32 


0 


40200 




FMU 


0 


THETADL 


40116 


20 


0 


40200 




STA 


0 


THETADL 




12 


0 


40664 




LDA 


0 


ONE 


40117 


30 


0 


40407 




FAD 


0 


RESTITUT 




32 


0 


40411 




FMU 


0 


RHO 


40120 


32 


0 


40214 




FMU 


0 


THETADM 




30 


0 


40200 




FAD 


0 


THETADL 


40121 


33 


0 


40414 




FDV 


0 


TOTINERT 




32 


0 


40412 




FMU 


0 


MOINERT 


40122 


20 


0 


40200 




STA 


0 


THETADL 




20 


0 


40176 




STA 


0 


V 


40123 


13 


0 


40176 




LAC 


0 


V 




32 


0 


40423 




FMU 


0 


B 


40124 


20 


0 


40175 




STA 


0 


VDOT 




13 


0 


40214 




LAC 


0 


THETADM 


40125 


32 


0 


40411 




FMU 


0 


RHO 




30 


0 


40212 




FAD 


0 


VJ 


40126 


32 


0 


40407 




FMU 


0 


RESTITUT 




30 


0 


40200 




FAD 


0 


THETADL 


40127 


33 


0 


40411 




FDV 


0 


RHO 




20 


0 


40214 




STA 


0 


THETADM 


40130 


20 


0 


40212 




STA 


0 


W 


" 


13 


0 


40212 




LAC 


0 


w 


40131 


32 


0 


40422 




FMU 
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PHASE TRAJECTORY PRINT OUT 

(Numbers indicated in tenths followed by power of 10 ) 
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